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ABSTRACT 

We analyse the results of a 5.5-yr photometric campaign that monitored 247 southern, semi- 
regular variables with relatively precise Hipparcos parallaxes to demonstrate an unambigu- 
ous detection of Red Giant Branch (RGB) pulsations in the solar neighbourhood. We show 
that Sequence A' contains a mixture of AGB and RGB stars, as indicated by a tempera- 
ture related shift at the TRGB. Large Magellanic Cloud (LMC) and Galactic sequences are 
compared in several ways to show that the P-L sequence zero-points have a negligible metal- 
licity dependence. We describe a new method to determine absolute magnitudes from pul- 
sation periods and calibrate the LMC distance modulus using Hipparcos parallaxes to find 
Mlmc = 18.54 ± 0.03 mag. Several sources of systematic error are discussed to explain dis- 
crepancies between the MACHO and OGLE sequences in the LMC. We derive a relative dis- 
tance modulus of the Small Magellanic Cloud (SMC) relative to the LMC of Afi = 0.41±0.02 
mag. A comparison of other pulsation properties, including period-amplitude and luminosity- 
amplitude relations, confirms that RGB pulsation properties are consistent and universal, in- 
dicating that the RGB sequences are suitable as high-precision distance indicators. The M 
giants with the shortest periods bridge the gap between G and K giant solar-like oscillations 
and M-giant pulsation, revealing a smooth continuity as we ascend the giant branch. 
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1 INTRODUCTION 

The existence of multiple period-luminosity (P-L) relationships 
is now well-established among Miras and semi-regular variables 
(SRVs) in the Large Magellanic Cloud (LMC), Small Magel- 
lanic Cloud (SMC), and Galactic Bulge. In the last decade, mi- 
crolensing surveys have provided a wealth of information, re- 
sulting in a co mpletely new picture of red giant pulsations 
(see, for example, Wood et al. 19991: iKiss & Bedd ing 2003 ,120041; 
Itaetal. 2004; Groenewegen 2004; Schultheis, Glass & Cioni 
20041 : IWrav. Ever & Paczvriskil2004l ; ISoszvnski et alj|2007l) . Most 
studies have concentrated on the (K, log P) plane, where M-giant 
emissions are strongest and interstellar extinction is relatively low, 
although similar P-L sequences have recently been shown for the 
mid-IR JGlass et all2009l) . 

M giants are common and intrinsically luminous, making 
them potentially important as distance indicators. Indeed, their P-L 
relations have already been used to measur e the distances to the ex- 
ternal galaxies NGC 5128 teeikubj|2004T) and the Phoenix Dwarf 
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JMenzies et alj 120080 . to determine the r elative distance modu- 
lus between the LMC, SMC, and Bulge fass & Bedding] I2004J: 
iGlass et al. 2009 ), and to probe the thr ee-dimensional structure of 
the LMC JLah. Kiss & Bed ding 2005). However, the metallicity 
dependence of the P-L relations is still uncertain, limiting their use- 
fulness. To determine whether the P-L relations are truly universal, 
we must first identify the local sequences in the solar neighbour- 
hood. 

The first evidence of m ultiple local sequences was provided 
by Bedding & Zijlstra J1998I) . who used Hipparcos parallaxes and 
published periods for 24 nearby M giants to de monstrate two dis- 
tinct seque nces. Subsequent inv estigations bv lKnapp et al.l (2003) 
and lYesilvaprak & Aslanl J2004I) confirmed at least two sequences, 
but these studies were limited by relatively large parallax errors, 
and a reduced sample s ize (by tightly c onstraining parallax errors), 
respectively. Recently, IGlass & van Le euwen ( 2007) used revised 
Hipparcos parallaxes and published periods for 64 red giants, con- 
firming that local stars obey similar P-L relations to those found 
in the Magellanic Clouds and Bulge, although details remained 
sketchy due to insufficient data, particularly for short-period stars. 

In Paper I JTabur et alj|20091) . we provided periods and K- 
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band magnitudes for 247 nearby red giants that we observed over 
5.5 years. With the benefit of revised Hipparcos parallaxes for this 
large sample, we now analyse those data to study the local se- 
quences. We describe our sample in Sect. [2] and provide period- 
luminosity diagrams in Sect. [3] The sequences are compared to 
those in the LMC to determine the effect of metallicity on the P-L 
zero-point (Sect.|4]&[5]l. A new method to estimate absolute magni- 
tudes from periods alone is described in Sect. [6] and is used to cal- 
ibrate the LMC distance modulus using Hipparcos parallaxes. We 
compare pulsation amplitudes to those in previous studies (Sect. 
|7]l and demonstrate how our sample relates to G and K giants that 
exhibit solar-like oscillations (Sect. [8}, and summarise our conclu- 
sions in Sect. [9] 



lower amplitudes mean the photometric errors become more signif- 
icant. 

Figure [2] shows normalised power spectra for the subset of 
stars with the most precise parallaxes (ov/V ^ 0.04), with Mk 
rounded to the nearest 0. 1 mag to avoid overlapping spectra. Dot- 
ted lines indicat e the positions of the P-L sequences using the 
nomenclature of lSoszvriski et al.l <2007. hereafter S07). These plots 
demonstrate a clear detection of the P-L sequences, and the preva- 
lence of short-period pulsation at low luminosities. 



2 DATA 

Our sample consists of 243 nearby M giantfjj with periods and 
absolute K magnitudes determined in Paper I, supp lemented by 
64 additional stars tabulated in Table 2 of iGlass & van Leeuwenl 
(2007). We used our periods in preference to the literature values 
for 21 stars that were in common. 

The periods derived from CCD observations in Paper I were 
limited to P < 300 d, since we considered those to be the most 
reliable. Longer periods were less well-determined due to gaps in 
the sampling, possible instrumental drifts, and difficulty match- 
ing light curves in different survey phases. However, periods de- 
rived using photoelectric p hotometry and those extracted from 
IGlass & van Leeuwen (2007) were retained for all values of the pe- 
riod. 

Apparent A'-band magnitudes were either determined from 
the best available NIR catalogue, as described in Pap e r I, o r 
obtained directly from Table 2 of IGlass & van Leeuwenl (2007). 
For HP 201298, the value of m Ks = 1.06 mag listed by 
IGlass & van Leeuwenl J2007I) is inconsistent wi th the value give n 
in the Two-Micron All-Sky Survey (2MASS: ICutri et alj|2003l) , 
namely rriK 3 = 1.729 ± 0.248 mag, and we hav e adopted the 
latter value. The use of revised Hipparcos parallaxes (van Leeuwenl 
120071) resulted in a reduction of the median relative parallax un- 
certainty of the sample from a-nj-n = 0.175 to 0.106, resulting in a 
median uncertainty in Mk of 0.23 mag. To sharpen the P-L ridges, 
we restricted our analysis to a subset of stars with the smallest rela- 
tive parallax uncertainties, as described in Sect.[3] However, before 
discussing the sequences of periods, we first present examples of 
the power spectra from which they were extracted. 



2.1 Power spectra 

Figure Q] shows a selection of Fourier spectra ordered by decreas- 
ing Mk- The ordinate shows the power, being the square of the 
Fourier semi-amplitude. The overall trend toward shorter periods 
and smaller amplitudes with decreasing luminosity is clear (note 
the change in vertical scale). Also visible is the correlation between 
adjacent spectra, where the dominant peaks are slowly shifting to- 
ward shorter periods as luminosity decreases and, in some cases, 
less prominent, shorter-period peaks start to dominate (see, for ex- 
ample, the transition from FH Eri to TT Crv, and NSV 18257 to 
NSV 15560). Noise increases toward the bottom of the plot, where 



1 We have excluded 4 stars with spectral types C or S from the original list 
of 247 stars with periods published in Paper I. 



3 PERIOD-LUMINOSITY RELATIONS 

Figure [3] shows period-luminosity diagrams in the (Mk, logP) 
plane for Galactic SRVs with five different thresholds on the par- 
allax uncertainties (triangles). For comparison, we have plotted 
the corresponding LMC sequences from M assive Compact Halo 
Object (MACHO) survey data analysed by iDerekas et al] (2006, 
hereafter D0 6), assuming a distance modulus of /ilmc = 18.50 
(Alves 2004). We show Optical Gravitational Lensing Experiment 
(OGLE) data from S07 for relative parallax errors larger than 10 
per cent. The sequences are labe lled A', A, B, C an d D, following 
the nomenclature established bv | Wpod et alj dl999T) . Sequence A', 
first identified by ISoszvhski et alj J20041) . consists of the shortest 
periods and smallest amplitudes, and is well-resolved in the sparser 
MACHO data. 

We briefly note a number of features in Figure [3] and discuss 
them in greater detail in following sections. Firstly, there is broad 
agreement between our Galactic sample and the LMC sequences, 
giving an overall impression of parallel ridges in both datasets. The 
agreement is excellent at small a w /n, particularly for Sequences 
A' and A, which are relatively narrow. The Galactic sequences be- 
comes less clear as parallax restrictions are relaxed, presumably 
due to relative parallax uncertainties, which contribute up to ~ 0.5 
mag of scatter. Error bars have been omitted for clarity, but we show 
the median uncertainty in each panel. 

Stars wit h periods longer than 300 d (log P ~ 2.5) were 
obtained from G lass & van Leeuwenl d2007h . A few have periods 
longer than those on Sequence D, and may be spurious. Periods 
that fall between Sequences C and D we att ribute to low-fre quency 
stellar noise, presumably from convection dKiss e t al. 2006). 

The significant decrease in stellar density at luminosities 
brighter than Mk ~ —6.8 marks the tip of the Red Giant Branch 
(TRGB), consistent with the recent finding that the local TRGB is 
detectable using revised Hi pparcos parallaxes and NIR p hotometry 
from 2MASS and DIRBE ( Tabur. Kiss & Beddindl2009l) . 

The Galactic sequences extend to similar low luminosities and 
short periods {Mk ~ — 4 and P ~ 10 d) as seen in the OGLE data 
for the LMC (bottom two panels of Figure [3J, with good correla- 
tion. Figure [4] shows a close-up of the lower-left region of the P-L 
diagram using MACHO and OGLE periods in the top and bottom 
panels, respectively. The sequence extending along the left edge to 
Mk ~ —4.8 and P ~ 6.7 d, traces out the limit on hig h-order 
radial pulsations imposed by the acoustic cut-off frequency (Wood 
2006). Sequences are shown using the nomenclature in S07. Blue 
triangles show all periods associated with 1 1 stars having at least 
one period significantly shorter than a4, which corresponds to Se- 
quence A'. These stars have large uncertainties in Mk, limiting the 
possibility of unique ridge identification. 
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Figure 1. A selection of power spectra ordered by decreasing absolute lumi- 
nosity that demonstrate the trend toward low-amplitude, short-period pulsa- 
tion with decreasing luminosity, giving the visual impression of sequences 
running diagonally down the plot. Significant periods (selected from Paper 
I) are marked with dotted vertical lines. Dashed horizontal lines indicate the 
mean signal in the frequency range 0.2-0.4 cycle d _1 . 
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Figure 2. Normalised power spectra of stars with <r n /n ^ 0.04 and Mk 
rounded to 0.1 mag. LMC sequences from S07 are shown for comparison, 
assuming /ilmc = 18.50. Dashed (dotted) lines represent Sequences 61— 
63 (al-a4) respectively. 



4 LMC SEQUENCES: MACHO VS OGLE 

To test for any metallicity dependence in the zero-point of the P-L 
relations, we wish to compare our data to the LMC, whose metal- 
licity is approximately half solar: Zq = 0.016 and Zlmc = 0.008 
dVassiliadis & Wooall994T) . As a first step, we compared the loca- 
tion of the LMC P-L (K a ) relations derived by S07 using OGLE-II 
and OGLE-III photometry, with the MACHO sequences from D06 
(Figure[4] top panel). 

To improve precision, S07 derived linear P-L relations using 
least-squares fits to density maps of the more numerous OGLE pe- 
riods which, they noted, produced slopes that were steeper than 
would have been obtained using linear fits. They identified four se- 
quences above the TRGB, labelled al-a4, that all extend below 
the TRGB, and three sequences, 61-63, that were only populated 
below the TRGB. The slight period sh ift that occurs at the TRGB 
JKiss & Bed ding 2003: llta et~alll20 04^ arises from the temperature 
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Figure 3. Period-Luminosity diagram of Galactic SRVs with ov/i" ^ 
0.025,0.05,0.10,0.15,0.25 (red triangles) and LMC stars (grey dots), 
with an assumed LMC distance modulus of 18.50. The top 3 panels use 
MACHO periods from D06. The bottom 2 panels use OGLE periods from 
S07, and show the overlapping luminosity range. 
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Figure 4. Top: P-L sequences for local SRVs with cttt/tt ^ 0.25 (red cir- 
cle) and LMC (MACHO) periods from D06 (black cross). Sequences from 
S07 are shown, together with 64 identified by this study. Bottom: Same as 
top panel, but using LMC(OGLE) periods. Stars with periods shorter than 
aA are indicated by blue triangles, including their counterparts on other 
longer-period sequences. These stars trace out the edge of the acoustic cut- 
off frequency. 



difference between first- and second-ascent red giants, and implies 
that stars below the TRGB consist mostly of RGB stars. Notably, 
S07 described Sequence aA as containing AGB stars only, base d 
on an earlier analysis of OGLE-II data bv lSoszvnski et al.l (2004), 
who found no evidence for a period shift at the TRGB. 

Sequence a4 corresponds to Sequence A' in D06 and, above 
the TRGB, we see good agreement in Figure [4] with the periods 
derived from MACHO data. However, the extension of aA below 
the TRGB does not match MACHO periods, which are shifted by 
<51ogP ~ -0.05 from aA. The MACHO data clearly show the 
expected shift at the TRGB, consistent in magnitude and direction 
with the shift between the longer-period sequences, implying that 
Sequence "64" does exist. Moreover, in Section[5] we confirm the 
presence of Sequence 64 in our sample of Galactic SRVs. 

An alternate representation of the LMC sequences is shown 
in Figure [5] where we compare the distribution of periods as a 
function of vertical offset from a reference sequence, chosen here 
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Figure 5. Vertical offsets of stars relative to 63, above (blue dotted) and be- 
low (red solid) the TRGB. The top panel shows LMC OGLE stars from S07, 
binned in 0.025 mag increments. All other panel are binned in 0.075 mag 
increments. The second panel shows LMC MACHO stars from D06. The 
remaining panels show the distributions for Galactic SRVs with a n /n < 
0.03, 0.06, and 0.09, respectively, assuming /ilmc = 18.50. The verti- 
cal dotted lines indicate the best-fit positions of Sequences 61-64 to LMC 
(MACHO) data, and the solid black line shows a continuous distribution 
function, as explained in the text. 



to be 63 (defined as m Ks = -3.72 log P + 17.66 by S07). The 
top panel shows OGLE periods from S07, selected with a variable 
S/N threshold, as explained in Section [6] and binned in units of 
0.025 mag. The second panel contains MACHO periods from D06, 
binned by 0.075 mag. The bottom 3 panels show distributions for 
our Galactic sample assuming /jlmc = 18.50, and are described 
in the next section. Peaks in the distributions indicate the positions 
of the sequences. Dotted vertical lines show the best-fit positions 
of Sequences 61—64 for MACHO data, which were determined by 
binning the data multiple times, using a variety of bins sizes and 
starting offsets. There is a small disparity between peak positions 
in the two datasets. Although Sequences 62 and 63 are clearly iden- 
tified in both datasets, OGLE peaks are shifted by about +0.1 mag 
and +0.04 mag, respectively, relative to MACHO. Both 61 and 64 
are ill-defined in the OGLE data, with relatively broad distribu- 
tions, but are clearer in MACHO. The centre of the 61 distribution 
is approximately aligned with MACHO, but 64 is offset by ~ +0.2 
mag. 

We investigated reasons for the discrepancy which, left un- 
solved, would limit possible conclusions about the effect of metal- 
licity. Several sources of systematic error were found to contribute. 

(i) Sample selection: The OGLE survey was restricted to the 
LMC bar, while MACHO observed a larger area, including regions 
outside the bar. To test for a selection effect, we obtained the origi- 
nal MACHO periodf|j corresponding to the subset of stars overlap- 
ping the OGLE sample. The 0.1 mag difference for Sequence 62 is 
no longer evident, with Sequences 61-63 closely aligned with their 
OGLE counterparts (Figure [6J. The sparsely populated sequence 
64 was not detected. Gaussian fits to the peaks show differences 
in their centres of 0.01 mag for 61 and 62, and 0.03 mag for 63. 
Although the MACHO sample contains fewer stars (its distribution 
has been multiplied by a factor of 50 in Figure[6j, the relative peak 
sizes approximate that of the larger OGLE sample. 

(ii) Systematics within the LMC bar: A small systematic error 
is introduced by sampling along the bar, which is inc lined ~ 30° to 
our line of sight, with the eastern end being closest dSubramaniuml 
l2003l : lLah. Kiss & Beddindl2005h . To quantify this, we divided the 
OGLE sample into 3 groups, each spanning 4° of RA, and calcu- 
lated distributions for each group (Figure |7). Gaussian fits to the 
peaks showed variations in relative distance modulus of 0.03-0.05 
mag. A sample consisting of a mixture of stars evenly distributed 
along the bar will tend to smooth out this effect. 

(iii) S/N thresholds: We selected samples from OGLE using sev- 
eral variable S/N thresholds, which improve the resolution of the 
sequences (as explained in Section [6), as well as a constant value 
of 7.5, used in Paper I. Sequence 63 showed a +0.03 mag sys- 
tematic shift as more stars were included, while 62, which contains 
fewer stars, showed a shift of +0.08 mag (Figure[8j. The location 
of 64 was very sensitive to the threshold. Using a variable thresh- 
old of 41ogP — 0.7, stars with periods on Sequence 64 had an 
effective S/N threshold of ~ 3, equivalent to the fixed S/N limit 
used by D06 to analyse MACHO data. This resulted in a shift of 
~ —0.15 mag compared to higher thresholds (which were unable 
to resolve the low-amplitude stars on this sequence), and explains 
the large difference between the MACHO and OGLE positions for 
this sequence. 

(iv) The frequency of stars occupying sequences may not be 
universal and may introduce a systematic error, particularly when 



We selected stars classified as LPVs from http://macho.anu.edu.au/star 
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Figure 6. The OGLE and MACHO sequences below the TRGB (thick black 
and thin red lines, respectively) are closely aligned when the MACHO sam- 
ple is restricted to the overlapping region of the LMC bar. The peaks are 
nearly coincident with the Galactic distribution (ov/ir < 0.03, blue dashed 
line), which is shown for comparison, assuming /ilmc = 18.50. The MA- 
CHO and Galactic distributions have been multiplied by factors of 50 and 
100, respectively, to show the relative sizes of the peaks. 
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Figure 7. Subsets of the OGLE data divided into 3 regions spanning 4° 
of right ascension along the LMC bar, binned in units of 0.025 mag. The 
peak offsets indicate a systematic error of ~ 0.03-0.05 mag occurs in the 
relative distance modulus, due to the tilt of the bar. 



using distribution-dependent analysis techniques, such as cross- 
correlation. In addition, choices made during processing, such as 
the S/N threshold, affects the distribution of stars within sequences. 
Figure [8] top panel, shows a large change in frequency of 64 rela- 
tive to 62 for different variable thresholds, while the bottom panel 
shows the effect of using fixed and variable cuts on Sequences 62 
and 63. We show that these moderate changes in distribution pro- 
duce relatively small systematic errors in Sect. [6] 

(v) Number of periods used: Restricting the number of period 
used for each star to a maximum of 1, 2, 3, 5 or 15 periods had a 
negligible effect on the positions of 62 and 63, which varied by less 
than 0.01 mag. 

We conclude that with appropriate sample selection and process- 
ing, the original MACHO and OGLE periods produce consistent 
results, and adopt the larger OGLE sample as being representative 
of the LMC sequences for the remainder of this analysis. 



5 EFFECT OF METALLICITY ON P-L ZERO POINT 

The slopes of LMC and Galactic P-L relations have been shown 
to be very similar (Fouque et al. 2007; Glass et al. 2009). Thus, a 
systematic difference in peak positions between the corresponding 
sequences would be interpreted as a difference in the zero-points 
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Figure 8. The distribution of stars in OGLE sequences for fixed and variable 
S/N thresholds. Top: A change in relative frequency and offset of 64 relative 
to 62 occurs when the threshold is low enough to include low-amplitude 
stars on 64. Bottom: A comparison of variable and fixed S/N cuts, showing 
the change in the number of stars in 62 relative to 63, due to the different 
responses to the rising noise floor at longer periods. 



of their P-L relations. To investigate the possibility, we show the 
Galactic M-giant distributions for stars with ov/V ^ 0.03, 0.06 
and 0.09, respectively (bottom three panels of Figure [5), To take 
into account the spread in the photometric and parallax uncertain- 
ties, we calculated a continuous distribution function for stars be- 
low the TRGB (solid black line) as a sum of Gaussians 
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where Oi and Oi are their offsets relative to 63 and their uncer- 
tainties, respectively, and N is the total number of periods in the 
sample JSakai et alJI 19960 . The agreement between the histograms 
and functions is good, with peaks closely aligned. Sequences 61, 
62, 63, and 64 (vertical dotted lines) are clearly identified at 3% 
and 6% relative parallax uncertainty, but less well-defined at 9%, 
due to blurring of the sequences. In light of our findings in the pre- 
vious section, we fitted Gaussians to the peaks of the Galactic se- 
quences (61, 62 and 63) and compared them to OGLE sequences 
defined by periods with S/N > 7.5, the same criterion used in 
Paper I. The difference in peak centres (OGLE — Galactic) ranged 
between -0.04 and +0.07 mag, and -0.08 and +0.09 mag for 
ov/tt ^ 0.03 and 0.06, respectively. The mean differences for the 
two samples were +0.03 ± 0.05 and +0.01 ± 0.07 mag, respec- 
tively, indicating identical peak centres to within the uncertainties. 
We used cross-correlation to assess the similarity of the OGLE 
and Galactic distributions below the TRGB, where the majority of 
periods were detected. A discrete cross-correlation function, 
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was applied, where x and y are the binned LMC and Galactic dis- 
tributions, respectively, N is the number of bins, and in is the 
trial offset. The best match is indicated by the value of m that 
maximises R. Since Figure [5] shows the distributions are closely 
aligned, only values of m that implied an offset less than 0.25 
mag were searched. R was calculated for a range of bin sizes and 
starting offsets to reduce sensitivity to these parameters, yielding 
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a mean offset and the rms scatter, which we adopted as our un- 
certainty. We found offsets implying that the Galactic distributions 
were fainter than their OGLE counterparts by 0.04 and 0.07 mag 
for samples with a w /n ^ 0.03 and 0.06, respectively, with uncer- 
tainties of ~ 0.03 mag. The expected K-band zero-point s hift for 
an in crease in metallicity by a factor of 2 is ~ 0.13 mag (Wood 
Il9901) . which is supported by observa tional evidenc e from a com- 
parison of SMC and LMC sequences Jlta et alj2004h . Although the 
Galactic shift is in the expected direction, its magnitude is signifi- 
cantly smaller than expected and is inconsistent with our previous 
result based upon peak centres. We caution that cross-correlation 
assumes an identical morphology of the signals being compared, 
which is not the case for the distribution of variables occupying the 
sequences in the two galaxies (Figure [6}. On the other hand, these 
comparisons assume ^lmc = 18.50 and, a small change in the 
adopted value would systematically alter the alignment of the dis- 
tributions, potentially affecting the conclusion. We use a novel new 
method in the following section to investigate further. 



6 ESTIMATING ABSOLUTE MAGNITUDES FROM 
PERIODS 

In order to use pulsating M giants as distance indicators, we now 
describe a new method that utilises the periods alone to determine 
the absolute magnitude of a given star, using the P-L sequences of 
the LMC as the reference. 

Paper I confirmed the multi-periodic nature of local M giants 
and highlighted the similarity of their pulsation properties to their 
counterparts in the LMC. Following the evidence presented in Sect. 
[5] we hypothesise that the LMC P-L sequences are representative 
of sequences in the SMC and local solar neighbourhood, with iden- 
tical slopes and zero-points. Using the distribution of LMC stars 
in the (mx s , log P) plane at specific periods, we derive probabil- 
ity density functions that describe the likelihood of a star having a 
particular magnitude. Using all measured periods for a given star, 
we estimate its magnitude by multiplying these probability density 
functions and locating the peak in the resulting probability func- 
tion. 

We constructed a fiducial LMC P-L diagram using OGLE- 
II and OGLE-III periods from S07. A^-band magnitudes were 
extracted from the Infrared Survey Facility catalogue (IRSF; 
iKato et al.N 2007). which provides deeper and more accurate pho- 
tometry than 2MASS, due in part to its superior resolution. To 
reduce the possibility of false matches, only stars that had a sin- 
gle match within a 1" radius were used. We transformed the IRSF 
magnitudes to the 2MASS system using the relation -Ks.irsf = 
Admass + 0.010(J - KT s ) 2 mass + 0.014 JKato et alj |2007h in 
order to compare the magnitudes, and found good agreement. How- 
ever, we noted a rapid divergence at faint magnitudes for stars with 
0.7 ^ log P < 0.9, which we excluded from the analysis. 

We combined OGLE data for OSARGs, Mira, SRV, and ellip- 
soidal variables into one dataset, while limiting the periods to the 
range 0.9 ^J log P ^ 3.5. Due to the rising noise-floor caused 
by 1// noise, we departed from the usual practise of using a fixed 
S/N cut-off. Instead, we used a sloping noise floor when calculating 
the S/N of each peak, rejecting periods with S/N < 5 log P — 0.7. 
This allowed only high S/N periods (at large P) to be selected, while 
still allowing short-periods with lower S/N to be used (limited by an 
additional constraint of S/N ^ 3). This scheme resulted in good 
resolution of all known sequences (Figure[9j. 

To estimate a star's magnitude, we first created a histogram of 
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Figure 9. LMC P-L diagram using OGLE periods and IRSF magnitudes. 
Periods were selected using an adaptive S/N cut-off which helped to resolve 
the shorter-period sequences at low S/N. Grey vertical bars correspond to 
the parts of the P-L used to create the probability density functions shown 
in Figure fTol 



the stellar density within a narrow, vertical strip of the P-L diagram 
(±2.5% wide in log P), centred about each of its periods. To pro- 
vide adequate resolution, counts were made for bins of 0.01 mag, 
which were then smoothed over ±5 bins using a boxcar average. 
Each resulting distribution, which is basically the luminosity func- 
tion of stars having that pulsation period, was converted to a prob- 
ability density function by normalising the area under the curve to 
unity. In this way, the function describes the probability of a star 
with that period having a particular magnitude. By combining mul- 
tiple periods by multiplying the corresponding probability density 
functions, an estimate of the magnitude of the star is given by the 
bin with the highest relative confidence. 

As an example, we consider a star in the LMC with 3 mea- 
sured periods: 46 d, 61 d, and 150 d. The corresponding vertical 
strips centred at log P = 1.66, 1.79 and 2.18 are shown in Figure 
[9] The normalised probability functions are shown in FigureQJJ to- 
gether with the combined probability, which indicates an estimated 
magnitude of m« s = 12.15. This is in close agreement with the 
measured value of m,K B — 12.18 ± 0.02. 

As a further test, we used the procedure to recreate the LMC 
P-L diagram using estimated magnitudes. The OGLE periods for 
each star used to produce Figure [9] were sorted by descending pul- 
sation amplitude. Magnitudes were estimated using the above tech- 
nique using the n most significant periods, where n = 1, 2, 3 and 
4 (Figure [Tit. A single period is sufficient to reproduce the main 
sequences, although the estimated mag at a given period always 
corresponds to the peak of the probability function and, thus, stars 
may be allocated to incorrect sequences. Adding further constraints 
(more periods) produces more reliable results, with the estimates 
no longer falling on the most populated sequence. Using 2 periods, 
Sequence A above the TRGB is resolved, and Sequences A' and B 
are more accurate. As the number of periods increases, a smooth- 
ing in the density of the P-L occurs, with gaps being filled in. We 
note that periods corresponding to faint magnitudes are moved up- 
ward in the diagram, due to the algorithm's propensity to gravitate 
toward areas of greater density in the P-L plane and, thus, it does 
not reproduce the known magnitudes in this case. This effect also 
explains the exaggerated gap at the TRGB. 

Figure [T2l shows the difference between the measured and es- 
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Figure 10. Probability density functions corresponding to the vertical strips 
in Figure [9] The combined probability (bottom right) indicates the most 
likely magnitude as irijt = 12.15, in close agreement with the measured 
value of 12.18 ± 0.02 mag. 
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Figure 11. LMC P-L diagrams using estimated magnitudes calculated us- 
ing \—A of the most significant periods for each star. The resolution of the 
sequences improves as more periods are used to estimate magnitudes. 



timated magnitudes, plotted against the confidence measure of each 
estimate, calculated as the height of the peak in the combined prob- 
ability distribution. The confidence is arbitrary, in the sense that it 
cannot be compared with other values produced with a different 
number of periods. It is only useful as a relative measure to com- 
pare estimates produced using the same number of periods. Figure 
[T2]confirms that magnitudes estimated from a single period have 
a significant chance of allocating stars to incorrect sequences, but 
results improve quickly as more periods are used to constrain the 
estimate. Estimates using 2, 3 and 4 periods produce progressively 
fewer outliers and a narrow distribution centred about zero, indicat- 
ing that the algorithm is working well. 
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Figure 12. LMC confidence measures as a function of magnitude difference 
(estimated — measured). The scatter decreases as more periods are used to 
estimate magnitudes. 



6.1 The SMC distance modulus 

We now apply this method to derive the distance modulus of the 
SMC relative to the LMC. As before, we obtained periods from 
OGLE-II and III (S07), and matched the SMC stars to the IRSF cat- 
alogue to obtain JiVband magnitudes. Magnitudes were estimated 
using the 2 or 3 strongest periods (largest amplitude), since sam- 
ple sizes were significantly reduced when requiring more periods 
satisfying our adaptive S/N cut-off. 

We used an iterative approach to account for the different lu- 
minosity ranges in the two galaxies. An initial pass using all peri- 
ods yielded an approximate relative distance modulus of A/i (a 0.4 
mag. Using this value, a second pass omitted periods that did not 
overlap the luminosity range of the LMC (specifically, those that 
were fainter than the 95th percentile of the distribution at a given 
period), resulting in an improved estimate and a smoother distribu- 
tion, with fewer outliers. The left-hand panels of Figure[T3]show the 
magnitude difference plotted against the confidence. Clearly, the 
difference is non-zero, allowing us to measure the distance modu- 
lus of the SMC relative to the LMC. 

The right-hand panels show the magnitude differences binned 
in units of 0.05 mag. It should be emphasised that the distribution 
is not a true Gaussian, since the wings contain some stars that were 
allocated to incorrect sequences. However, the core of the peak is 
composed of stars with high confidence levels that were likely to 
be assigned to the correct sequence. Gaussian fits to the distribu- 
tions yield a relative SMC distance modulus of A/i = 0.39 ± 0.01 
mag and 0.40 ± 0.01 mag for 2 and 3 periods respectively, which 
are consistent with the values derived from the subset of stars with 
high confidence levels only. The quoted uncertainties are the formal 
uncertainties from the Gaussian fit. 

To assess the sensitivity of this method to systematic errors 
introduced by different S/N thresholds and, indirectly, the rela- 
tive number of periods in each sequence, the samples were repro- 
cessed using periods selected with the same fixed and variable S/N 
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Figure 13. Magnitude differences plotted against confidence (left panels) 
for SMC stars using 2 and 3 OGLE periods, respectively, and the corre- 
sponding distributions (right panels). The offset of ~ 0.40 mag represents 
the SMC distance modulus relative to the LMC. 



thresholds shown in Figure [8] The results were very similar, with 
a mean relative distance modulus of A/i = 0.40 ± 0.01 mag and 
0.41 ± 0.02 mag for 2 and 3 periods, respectively. Using an ex- 
tensive list of distances published since 2 002 to calculate a mean 
distance modulus for the LMC and SMC, iTammann et alj (2008) 
found a difference of 0.40 ± 0.02 mag, in excellent agreement with 
our result. 



6.2 Calibration of the LMC distance modulus using 
Hipparcos parallaxes 

We now apply the algorithm to our sample of nearby Galactic M 
giants to give the magnitude of each star relative to the LMC fidu- 
cial. Paper I listed all periods derived from our Fourier analysis, 
which included some closely spaced frequencies that were prob- 
ably caused by amplitude and phase modulations of a single os- 
cillation mode. To remove them, all periods for a given star were 
sorted by descending amplitude. Each period was then compared 
to the remaining periods (of lower amplitude) to remove those with 
a period ratio < 1.1, which were assumed to be phase/amplitude 
modulations of the primary period. Periods longer than 150 d were 
also excluded, to avoid spurious periods caused by 1// noise. 

To improve the statistics, we used all stars in the sample hav- 
ing 2-5 periods, regardless of their relative parallax and photomet- 
ric uncertainties. We propagated uncertainties in the observed pe- 
riods, assumed to be 3%, to estimated magnitudes using a Monte- 
Carlo simulation with 10000 iterations for each star. In most cases, 
the estimated magnitudes derived from the observed periods were 
close to the simulation values (Figure[l4] top left). Stars with large 
uncertainties were found to have either broad peaks or multiple 
peaks of similar height in their probability function, which yielded 
disparate estimates from only small changes to the input periods. 
We show the combined uncertainties from period, parallax and pho- 
tometric uncertainties in Figure [T4j top right. 
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Figure 14. Galactic M giant results using 2-5 periods. Top left: differ- 
ence between magnitude estimates derived from the observed periods and 
a Monte-Carlo simulation. Magnitude differences (Mx - estimated) as a 
function of uncertainty (top right) and confidence (bottom left). Bottom 
right: distribution of estimated magnitudes for the full sample (dashed) and 
the subset with uncertainties < 0.15 mag (dotted). 



We used two passes to determine the LMC distance modu- 
lus. Due to the use of a variable number of periods per star, con- 
fidence values are not easily comparable, although a logarithmic 
plot shows four distinct groups, indicating the number periods used 
(Figure[l4] bottom left). The distribution of magnitude offsets with 
respect to the LMC is shown in the bottom right panel, binned by 
0.3 mag (dashed line). The distribution for the subset of stars with 
uncertainties smaller than 0.15 mag is shown with a dotted line. 
Gaussian fits to the two distributions (solid lines), yield an identi- 
cal distance modulus of /zlmc = 18.54 mag, with a formal un- 
certainty of ~ 0.02 mag. Reprocessing our Galactic sample with 
OGLE periods selected with the S/N thresholds shown in Figure [8] 
yielded a relatively small range of distance modulus spanning 0.08 
mag, with a mean of /ilmc = 18.54 ± 0.03 mag. Tests using the 
subset of stars with the smallest uncertainties produced nearly iden- 
tical results. Our distance modulus agrees with the mean value of 

18.53 ± 0.01 mag derived from 16 r ecent, independent, distance 
determinations dTammann et alj |2008), which is identical to within 
the uncertainties. Thus, our data show that the P-L relation zero- 
points have a negligible metallicity dependence. 

As a final test of this method, we reprocessed our sample with 
OGLE periods selected from stars within the three R.A. ranges 
shown in Figure[7] in an attempt to detect a systematic shift caused 
by the tilt of the LMC bar. This yielded values of ^ilmc = 18.59, 

18.54 and 18.52 ± 0.01 mag for ranges centred on R.A. = 77°, 
81° and 85°, respectively. These distances correctly show that the 
eastern end of the bar is closest, while the range in relative distance 
modul us agrees with th e recently determined value of 0.1 ± 0.03 
mag dLah. Kiss & Bedd ing 2005). 

Our method has some advantages over the Mira P-L relation- 
ship for determining distances. SRVs are more numerous and, their 
multi-periodic nature places constraints on their absolute magni- 
tude, allowing a single star to be used to estimate distance. Miras 
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are more luminous and have higher amplitudes, placing less de- 
mands on photometric precision. Both methods require monitoring 
over similar time scales to measure either one long period or, multi- 
ple shorter periods. However, unlike short-period SRVs, Miras are 
more affected by circumstellar extinction due to dust, which ca n 
be significant at periods longer than ~ 125 d JAlard et alj|200ll) . 
and introduces an additional source of uncertainty. Our method 
makes use of both Mira and SRV sequences and, having a negli- 
gible metallicity dependence, is well-suited to measuring distances 
in the solar neighbourhood, Galaxy and beyond. 



7 PULSATION AMPLITUDES 

The usefulness of M giants as distance indicators relies upon RGB 
pulsation properties being universal. Paper I showed that both pe- 
riod ratios and the incidence of multi-periodicity in nearby M gi- 
ants are consistent with stars in the Galactic Bulge and Magellanic 
Clouds. We now examine the pulsation amplitudes of nearby M gi- 
ants, to better understand the limitations, if any, they impose on our 
understanding of the P-L sequences. It is also interesting to study 
pulsation amplitudes in their own right, for example, to determine 
the relative roles of Mira-like and solar-like excitation processes 
(see Sect. [8). 

Our survey detected peak-to-peak Fourier amplitudes ranging 
from 0.86 mag (ir 1 Gru) down to 8 mmag (<f> 2 Hya), with 90% 
of periods having amplitudes A < 0.1 mag, and 75% having A < 
0.05 mag. We show period-amplitude and luminosity-amplitude di- 
agrams for stars with ov/V ^ 0.25 in Figure Q3] Different sym- 
bols have been used for periods corresponding to fundamental- 
mode pulsation (labelled C), higher overtones (labelled A,B), long 
secondary periods which have unknow n origin; see, for example, 



IWood et al.ll2004 iNicholls et alj|2009l (D), and those falling be- 
tween Sequences C and D (labelled C-D). Periods were assigned 
to particular sequence when they were with in ±0.5 mag of the re- 
lations defined by S07 (Table 1). Data from lGlass & van Leeuwenl 
(2007) were omitted since amplitude data were not available for all 
published periods. 

In the period-amplitude diagram (left panel), the majority of 
detected periods correspond to short-period overtone pulsation. A 
clear trend of increasing pulsation ampli tude with period is evident, 
including a fairly sharp upper envelope. lAlard et alj ( 1200 lb studied 
SRVs in the low-extinction region of Baade's window of the Galac- 
tic Bulge, finding that amplitude increased with period, reaching a 
typical value of about 0.3 mag at 100 d, in agree ment with our re- 
sults. Similarly. IWrav. Ever & Paczvfiskil J2004l) identified 15000 
red variables in the Galactic bar, with a typical /-band amplitude 
of about 0.2 mag at 100 d. An increase in pulsation amplitude with 
lumino sity, and therefore in period , is expected for solar-like oscil- 
lations JKieldsen & Beddindll995l) . The observations for overtone 
pulsators (A and B) lend support to the idea that pulsations in these 
stars are excited stochastically by convection. 

Sequence C shows a weaker correlation between period and 
amplitude, in agreement with the results of D06 for the LMC. Se- 
quence D contains very few periods, but a positive correlation with 
amplitude is also visible. Notably, the periods falling between Se- 
quences C and D, which we have attributed to low-frequency noise 
(see Sect. [2j, are positively correlated with amplitude and have a 
slope similar to the other sequences. This result suggests a physical 
origin and can be explained by stellar noise from convection, which 
is expected to increase towards lower frequencies. OGLE detected 
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Figure IS. Period-amplitude (left) and luminosity-amplitude (right) rela- 
tions for our survey fields. CCD amplitudes are shown in the instrumental 
system, which has a broad response that peaks betwe e n V a nd R bands. 
Squares denote periods on Sequence C of Wood et al. (1999), which con- 
tains stars pulsating in the fundamental mode, while crosses are used for 
those pulsating in overtones. Open circles denote stars on seq D, contain- 
ing stars with Long Secondary Periods, and triangles denote stars between 
Sequences C & D. 



a significant number of periods between Sequences C and D (see 
Figure|9](, consistent with our result. 

The luminosity-amplitude diagram (right panel of Figure IT3T> 
shows that pulsation amplitudes are positively correlated with ab- 
solute luminosity for all sequences. A sudden decrease in density 
occurs arou nd Mk ~ —6.85 mag ( dashed line), which marks 
the TRGB JTabur. Kiss & Beddindl2009l) . With the exception of 
a few outliers, the majority of pulsation-related amplitudes below 
the TRGB are smaller than 0.2 mag, in agreement with the period- 
luminosity-amplitude data-cube slices of the LMC and SMC, 
which show very few stars with A > 0.14 mag jKiss & Bedding 
I2003ll2004l) . 



7.1 V1472 Aquilae 

Stars in the luminosity-amplitude plot obey a well-defined ampli- 
tude limit at a given luminosity, with the exception of V 1472 Aql, 
whose amplitude is unexpectedly large for its absolute magnitude 
(M K = -3.78, A=0A2 mag, P=100.03 d), suggesting that this pe- 
riod is not pulsation related. Our period and amplitude determina- 
tions are in close agreement with the Hipparcos values of P=100.37 
d and ^4h p =0.16 mag, both of which place the star on the low- 
luminosity exten sion of Sequenc e D, w hich contains ellipsoidal 
variables (D06). JLucke & Mavorl dl982l) used 41 radial velocity 
measurements to determine a spectroscopic period of 198 d, and 
concluded that it could be a contact binary that may show eclipses. 
Indeed, using the spectroscopic period, ISamusI ( 119970 folded the 
Hipparcos photometry to plot a low-quality light curve revealing 
continuous brightness variations, and possible indications of pri- 
mary and secondary minima, consistent with an eclipsing or ellip- 
soidal variable star. Our data show that V1472 Aql exhibits the pri- 
mary and secondary minima of an ellipsoidal system (Figure [7(51). 
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Figure 16. Phased light curve of V1472 Aql folded with P=200.05 d. The 
continuous brightness variations and primary and secondary minima are 
typical of an ellipsoidal variable. 



8 BRIDGING THE GAP TO G AND K GIANTS 

The most luminous AGB stars undergo stable, mono-periodic 
Mira-like pulsations that are driven by opacity changes. On the 
other hand, solar-like oscillations, the stochastically driven ex- 
citation and damping of linearly stable modes by convection 
dDziembowski et al. 2001) may explain the less regular light vari - 
ations of SRVs dChristensen-Dalsgaard. Kieldsen & Matteill200ll) . 
Examples of sol ar-like oscillatio ns have been identified amongst 
Galactic SRVs dBeddingJ 120031) . including L 2 Pup, which ex- 
hibits characteri stics of both Mi ra-like pulsation and solar-like 
oscillations ( [Bedding et alj 12005ft . Low-amplitude, short-period 
photometric variability has also been observed in G and K gi- 
ants, with amplitudes of 5-15 mma g and periods of 2-4 d 
fedmonds & Gill iland 1996: lHenrv~e7a l. 2000), leading to the con- 
clusion that pulsation may extend to spectral classes earlier than 
M. However the precise nature of the driving mechanism is still 
unknown. 

Recently, solar-like oscillations have been observed in a num- 
ber of G and K giants with well determined parallaxes, allowing 
them to be added to our P-L diagram. Table Q] summarises the re- 
sults from 7 recent papers, tabulating the Hipparcos catalogue num- 
ber, star name, spectral type, revised Hipparcos parallax and un- 
certainty ( Ivan Leeuwenll2007l) . apparent K magnitude, absolute K 
magnitude, magnitude source (2=2MASS, D=DIRBE, G=Gezari, 
m=mean of all three) as discussed in Paper I, period corresponding 
to the frequency of maximum power, and a reference to the original 
paper. Additionally, data for UK g iants observed wit h the WIRE 
satellite were taken from Table 1 of lStello et alj (12008). Periods in 
hundreds more G and K giants have been measure d by the space 
missi ons CoRoT JHekker et alj 12009!) and Kepler dBedding et alj 
120101) . However, with the exception of HR 7349, a lack of paral- 
laxes prevents them from being included in Figure [TTI 

The combined P-L diagram (Figure [T7J shows the LMC se- 
quences from S07 (grey circles), our local SRVs with ov/V ^ 0.25 
(red circles) and the G and K giants (blue triangles). The M giants 
have sufficiently short periods and low luminosities to bridge the 
gap to the most luminous, longest period G and K giants. 

The left-most edge of the M-giant sequences traces the acous- 
tic cutoff frequency, whose approximate location has been plotted 
by scaling the Sun's frequency of f ac ~ 5.5 mHz using the relation 
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(Bro wn et alJI 19911) . using stellar 
parameters fromlCarroll & Ostlie (2007) and V — K colours from 



lBessell&BretiUl988h . 

The emerging picture is a clear progression of properties as 



Table 1. G/K giants with relatively precise parallaxes that exhibit solar-like 
os cillations. Reference s: A-Frandsen et al. 12002), B=Retter et al. (2003), 
C ^Tarrant et all J2007J) . D JTarrant et alj J2008|) . E ^De Ridder et alj (2006). 
F JBarban et all 12004) . G ^arrier etaTl 120101) . 
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Figure 17. P-L diagram extended to include nearby Galactic G and K gi- 
ants exhibiting solar-like oscillations. A dashed line marks the approximate 
location of the acoustic cutoff frequency, scaled from the solar value of 5.5 
mHz, as explained in the text. 



we ascend the giant branch. Stochastically excited, low-amplitude, 
micro- variability in G and K giants progresses to a mixture of Mira- 
like pulsation and solar-like oscillations in SRVs, and ends with 
stable, large-amplitude, Mira-like pulsation. Our data nicely bridge 
the gap between G/K and M giant pulsation, demonstrating that a 
smooth continuity exists. 



9 CONCLUSION 

We have conducted a 5.5-yr survey of 243 southern, pulsating M 
giants in the local solar neighbourhood. Combining our data with 
additional literature periods, we have analysed a sample of 286 stars 
with relatively precise Hipparcos parallaxes to unambiguously de- 
tect 4 P-L sequences below the TRGB. Together with Paper I, we 
have shown that the pulsation-related properties of nearby M gi- 
ants, including period ratios, amplitudes, multi-periodicity, and the 
temperature related shift at the TRGB are consistent with SRVs 
found in the Magellanic Clouds and Bulge, indicating that pulsation 
in a truly universal and consistent phenomenon. The metallicity de- 
pendence of the P-L zero points are negligible over the metallicity 
range spanned by the SMC, LMC and Galaxy, indicating that SRV 
P-L relationships are suitable as high-quality distance indicators. 
A new method for estimating absolute magnitudes from pulsation 
periods is described, which yields an LMC distance modulus of 
Mlmc = 18.54 ±0.03 mag, calibrated using Hipparcos parallaxes. 
We demonstrate that V1472 Aql is an ellipsoidal variable, whose 
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period and absolute magnitude place it on the low-luminosity ex- 
tension of Sequence D, in agreement with earlier conclusions about 
the nature and corrected location of Sequence D. We show that 17 
nearby G and K giants that exhibit solar-like oscillations have peri- 
ods similarly constrained by the acoustic cutoff frequency. Our data 
bridge the gap, both in luminosity and period, between G and K gi- 
ant solar-like oscillations, and M giant pulsation, demonstrating a 
clear continuity as we ascend the giant branch. 
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